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Abstract 



\ In this paper, we introduce the notion of discrete embedding which is an algebraic 

procedure associating a numerical scheme to a given ordinary differential equation. We 
first define the Gauss finite differences embedding. In this setting, we study variational 
integrator on classical Lagrangian systems. Finally, we extend these constructions to 
the fractional case. In particular, we define the Gauss Griinwald-Letnikov embedding 
and the corresponding variational integrator on fractional Lagrangian systems. 
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in 

O ■ Introduction 

cn 

Fractional calculus is a source of many investigations for the purposes of Physics and 
Engineering (see [1], |15j . [18], |22) . |24j). For example, fractional terms appear numerically 
in the modeling of the flow of a fluid in a heterogeneous environment, |25| . In Mathematics, 
we find fractional calculus in the probability theory, |20| . Moreover, due to their nonlocal 
characteristic and consequently due to their capacity of memory, fractional differential 
operators are also interesting for financial mathematics, |10) . and for modeling of some 
materials like gum and rubber. 

In recent years, many studies have been devoted to fractional Lagrangian systems, |llj . 
They arise for example in fractional optimal control theory ([I], [2], |14j . |16j ) and they are 
difficult to solve explicitly (see |22| for a general study). Then, it is interesting to develop 
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efficient numerical schemes to such systems. 

There exists a suitable method for classical Lagrangian systems called variational in- 
tegrators which is developed in |17) . |21] , Indeed, classical Lagrangian systems possess a 
variational structure, i.e. their solutions correspond to critical points of a functional and 
this characterization does not depend on the coordinates system, [5]. This structure in- 
duces strong constraints on the qualitative behavior of the solutions, as for example the 
conservation of energy for autonomous classical Lagrangian systems. 

The basic idea of a variational integrator is to preserve this variational structure at the 
discrete level. We can obtain this numerical scheme by forming a discrete analogous of the 
variational principle on a discrete version of the Lagrangian functional. 

Following the usual strategy of |17j and |21| . this paper is devoted to the extension of 
variational integrators to the fractional case in the framework of discrete embeddings. 

More precisely, we first define the notion of discrete embedding which gives a direct 
discrete analogous of a given differential equation (in particular of a given Lagrangian 
system) . This procedure is algebraic and mainly based on the form of differential operators 
which is coordinates dependent. Consequently, this process does not necessary conserve at 
the discrete level the intrinsic Lagrangian structure of a Lagrangian system. 

On the other hand, a discrete embedding gives a discrete version of a given Lagrangian 
functional and we can develop a discrete calculus of variations on this one which leads to 
a variational integrator of the associated Lagrangian system. 

Then, by defining a discrete embedding, we obtain two discretizations of a Lagrangian 
system: the direct one and the variational integrator. Of course, these two discretizations 
are not necessary the same: in this case, we say that the discrete embedding is not coher- 
ent. 

In this paper, we first prove that the basic finite differences discrete embedding is not 
coherent in the classical case. However, this default of coherence can be corrected by 
rewriting, with left and right derivatives, the classical differential equation of Lagrangian 
systems. Indeed, this asymmetric rewriting corresponds at the continuous level to the dis- 
crete version of Lagrangian systems obtained by variational integrator. As a consequence, 
the discrete embedding is coherent in the asymmetric case. 

In the fractional case, the derivatives are naturally asymmetric and satisfy a fractional 
integration by parts. In this paper, we prove that the fractional Griinwald-Letnikov dis- 
crete embedding is coherent. 

The discretization of fractional Euler-Lagrange equations has been discussed by several 
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authors. We refer in particular to [3] for finite element methods to fractional Lagrangian 
functionals, to [7] and [5] for a discrete fractional calculus of variations and to [6] for a 
numerical scheme on fractional optimal control problems. Some preliminary results on 
fractional discrete operators have already been discussed in these papers. 

However, the discrete embedding point of view and the associated notion of variational 
integrator are not introduced in these papers as well as the corresponding coherence prob- 
lem. 



The paper is organized as follows. In section [T] we define the notion of discrete embed- 
dings and direct discrete embeddings of a differential equation. Section [2] recalls definitions 
and results concerning Lagrangian systems. Then, we apply the previous theory of discrete 
embeddings to Lagrangian systems. Section [3] recalls the strategy of variational integrators 
in the framework of discrete embeddings. Noticing a default of coherence of the previous 
discrete embedding, we define in section U] asymmetric derivatives and then asymmetric 
Lagrangian systems. With these asymmetric notations, we correct the default of coher- 
ence. Section [5] is an introduction to fractional calculus. Finally, section [6] is devoted to 
fractional discrete embeddings and the associated fractional variational integrators. 



1 Notion of discrete embeddings 

In this paper, we consider classical and fractional differential systems in M. d where d 6 N* is 

the dimension. The trajectories of these systems are curves q in C°([a, 6],R a! ) where a < b 

dq d^q 

are two reals. For smooth enough functions q, we denote q = — and q = — 

dt dt 1 

1.1 Discrete embeddings 

A discrete embedding is a particular way to give a discrete analogous of an ordinary dif- 
ferential equation: 

Definition 1. Defining a discrete embedding means giving a discrete version of the fol- 
lowing elements: the curves q G C°([a, fe],R ), the derivative operator — and the functionals 
a : C°([a, b],M. d ) — > R. More precisely, it means giving: 

• an application q i — > q h where q h G (M. d ) mi , 

• a discrete operator A : (IR d ) mi — y (M. d ) m2 which discretizes the differential operator 
d 

dt' 

• an application a i — > a h where a h : (W 1 )™ 11 — > R, 
where mi , m2 € N* . 
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Obviously, the discretization of curves above-defined induces the discrete evaluation of a 
function / € C°(R d ,M m ), m € N* on the curves q. Indeed, by abusing the notation, if we 
denote by / the following application 

/: C°([a,b],R d ) — > C°([a,b],R m ) , 
q i — > f oq 

then the discretization of / is f h given by: 

f h : (R d ) mi — > (R m ) mi 

q h _> f h {qh) := (/(«£)) 

\ / fc=l,...,mi 

However, the discretization of a functional is not obvious. Indeed, it could use the dis- 
cretization of — and the discretization of other mathematical tools depending on the form 
dt 

of the functional itself. We are giving an example in subsection 13.11 with the discretization 
of a Lagrangian functional is written with an integral. 

In order to illustrate definition [TJ we define forward and backward finite differences em- 
bedding. For all the rest of the paper, we fix a = ± and N € N*. We denote by h = ^ 
the step size of the discretization and r = (tk)k=Q ... N the following partition of [a, b\. 

Vk G {0, N}, t k = a + k b -^. 
Now, we can define the following discrete embeddings: 

Definition 2 (case a = +). We call forward finite differences embedding denoted by 
FDE+ the definition of the following elements: the application 

disc: C°{[a,b],R d ) — > (R d ) N+1 

q i — > {q{t k )) k =o,...,N 

and the discrete operator 

A+: (R d f +1 — > (R d ) N 

r\ _ ( r\ \ , I Qk-Qk+i 

W — Wk)k=0,...,N 1 > ( h 



h J k=0,..,N-l 

Definition 3 (case a = — ). We call backward finite differences embedding denoted 
by FDE— the definition of the following elements: the application 

disc: C°{[a,b],R d ) — ► (R d ) N+1 

q i — ► (q(tk))k=o,...,N 

and the discrete operator 

A_: (R^)^ 1 — > {R d ) N 



Q = (Q 



k)k=0,...,N 



k=l,..,N 
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We use these notations in order to be uniform with fractional notations (in sections [MS]) . 

Let us notice that the discrete analogous of — in FDEa is then — o~A a . 

6 dt 

1.2 Direct discrete embedding 

Defining a discrete embedding allows us to define a direct discrete version of a differential 
equation: 

Definition 4. Let be fixed a discrete embedding defined as in definition^ and let (jEj) be 
an ordinary differential equation of unknown q E C°([a,b],M. d ) given by: 

O(q) = (E) 
d i 

where O is a differential operator shaped as O = ^J/i(0^ ° 9i(-) w here fi, gi are 

i 

functions. 

Then, the direct discrete embedding of (JEj) is (Eh) the system of equations of unknown 
q h € (M d ) mi given by: 

O h {q h ) = (E h ) 

where O h is the discretized operator of O given by O h = //*(•) A* ogf(.). 

i 

Example 1. We consider the Newton's equation with friction of unknown q G C°([o, b], M. d ) 
given by: 

Vt € [a,b], q(t) + q(t) + q(t) = 0. (NE) 

Then, the direct discrete embedding of ([NE]) with respect to FDE+ is ( NE^ ) the system 
of equations of unknown Q € (M <i ) Ar+1 given by: 

wi r- rn at oi Qk+2 - 2(Jfc+i + Qk Qk+i-Qk ^ n , ATTr 
vk G {0, N — 2}, — 2 1 \-Q k = 0. (NE h ) 

The direct discrete embedding of an ordinary differential equation is strongly depen- 
dent on the form of the differential operator O (and not on its equivalence class). The 
process O — > O h is not an application. For example, the discretized operator O h of 
O = — o sin{.) = cos(.) is different depending on the writing of O. 

\JjL (Jib 



2 Discrete embeddings of Lagrangian systems 
2.1 Lagrangian systems 

In this section, we recall classical definitions and theorems concerning Lagrangian systems. 
We refer to [5 J for a detailed study and for the proof of theorem [TJ 



Definition 5. A Lagrangian functional is an application defined by: 

L: C 2 ([a,b},M d ) — ► R (2.1) 
q i— ► / L(q(t),q(t),t)dt 

J a 

where L is a Lagrangian i.e. a C 2 application defined by: 

L : R d x R d x [a, b] — ► R 

(x,u,i) i — )• L(x,v,t). 

Let £ be a Lagrangian functional, we denote by DC(q)(w) the Frechet derivative of C in 
q along the direction w in C 2 ([a,b],M. d ), i.e. 

g/:( g )H = iim £(g + £w) - £(g) . 

e-K) £ 

An extremal (or critical point) of a Lagrangian functional £ is a trajectory q such that 
DC(q)(w) = for any variations w (i.e. u; G C 2 ([a, 6], R d ), u>(a) = u;(6) = 0). 
Extremals of a Lagrangian functional can be characterized as solution of a differential 
equation of order 2 given by: 

Theorem 1 (Variational principle). Let C be a Lagrangian functional associated to the 
Lagrangian L and let q G C 2 ([a,b],M. d ). Then, q is an extremal of £ if and only if q is 
solution of the Euler- Lagrange equation given by: 

VtE[a,b], ^(q(t),q(t),t) - ^(|(g(t),?(t),i)) =0. (EL) 

2.2 Direct discrete embedding of the Euler-Lagrange equation 

In this subsection, we apply definitions of section [T] in order to define the direct discrete 
embedded Euler-Lagrange equation. 

Definition 6. Let L be a Lagrangian and jEL]) its associated Euler-Lagrange equation. 
The direct discrete embedding of (lELp with respect to FDEa is given by: 

— (Q, -aA a Q, t) + aA a ( — (Q, -aA a Q, r) \ = 0, Q G (R d ) N+1 . (2.2) 

We illustrate this result with the classical example of the mechanical Lagrangian: 

L(x,v,t) = ^v 2 - U(x), (x, v, t) G R d x R d x [a, b] (*) 

where U represents the potential energy of the system. The mechanical Lagrangian gives 
the following Euler-Lagrange equation: 

Vt€[o,b], q(t) = -VU(q(t)), q€C 2 ([a,b],R d ). 
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Then, by direct discrete embedding with respect to FDE-, we obtain the following 
numerical scheme : 

\/ke2,...,N, Qk ~ 2Qk - 1 + Qk - 2 = -VU(Q k ), Q€(R d ) N+1 - 

3 Discrete embeddings and variational integrators on Lagrangian 
systems 

As we said previously, a direct discrete embedding is only based on the form of the dif- 
ferential operator which is dependent of the coordinates system. Then, arises the natural 
question of conservation of intrinsic properties of the differential equation at the discrete 
level. This paper is devoted to the conservation of the Lagrangian structure. In this section, 
we introduce variational integrators in the framework of the discrete embedding FDEo. 

3.1 Discrete calculus of variations and discrete embeddings 

In this subsection, we recall the strategy of variational integrators which are discretizations 
of Lagrangian systems preserving their Lagrangian structures. It consists in building the 
discrete analogous of the variational principle on a discretized Lagrangian functional. For 
more details, we refer to |17] and |21j . In our case, the discrete Lagrangian functional is 
obtained by giving a discrete embedding as defined in section [TJ 

3.1.1 Discrete Lagrangian functionals 

Giving FDEa induces the discretization of a Lagrangian functional as long as a quadrature 
formula is fixed in order to approximate integrals (as in (|2.ip ). We choose the usual o- 
quadrature formula of Gauss: 

for a continuous function / on [a, b], we discretize / f(t)dt by 

J a 

N-l N-l 



k=0 k=0 

and: 

N N 

Y,(tk - tk-i)f(t k ) =hJ2 f(tk) if o- = -. 

k=l k=l 

This process defines the Gauss finite differences embedding which we denote by 
Gauss-FDEa. We can now define the following discrete Lagrangian functional: 

Definition 7. Let C be a Lagrangian functional associated to the Lagrangian L. The 
discrete Lagrangian functional associated to C with respect to the discrete embedding 
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Gauss-FDEa is given by: 

C u h : (R d ) N+1 — ► R 

Q = (Qk)k=o,...,N 1 — > h ^ L(Qk, (—aA a Q)k,tk) 

k£l a 

where I + = {0, N - 1} and J_ ={1,...,JV}. 

Example 2. Lei us consider the Lagrangian functional associated to the mechanical La- 
grangian given by: 



C: C 2 ([a,b] 



, \ bl -q{tf-U{q{t))dt. 



Then, the discrete Lagrangian functional associated to C with respect to the discrete em- 
bedding Gauss-FDE— is given by: 

,N+1 



C h : (R d )^ — > R (3.1) 

Q = (Qk)k=o,...,N — ► hp^ ^-^ y -U(Q k ). 

3.1.2 Discrete calculus of variations 

Let £ be a Lagrangian functional and the discrete Lagrangian functional associated 
with respect to Gauss-FDEa. A discrete extremal (or discrete critical point) of 
is an element Q in (R^)-^" 1 " 1 such that DCl (Q)(W) = for any discrete variations W 
(i.e. W € (R d ) N+1 , W = W N = 0). 

Discrete extremals of can be characterized as solution of a system of equations: 

Theorem 2 (Discrete variational principle). Let be the discrete Lagrangian functional 
associated to the Lagrangian L with respect to Gauss-FDEa. Then, Q in (R d ) N+1 is a 
discrete extremal of £^ if and only if Q is solution of the following system of equations 
(called discrete Euler- Lagrange equation) given by: 

— (Q, -aA a Q, t) - aA_ CT I — (Q, -aA a Q, r) = 0, Q e (R d ) N+l . (EL h ) 

Thus, the variational integrator with respect to Gauss-FDEa on ([EL]) is the numerical 



scheme flEX/J : it corresponds to the discrete variational principle on the discrete La- 
grangian functional defined with respect to Gauss-FDEa. 

Example 3. As in example^ let us consider the Lagrangian Q. The Gauss-FDE— 
leads to the discrete Lagrangian functional given by (|3.ip . The discrete variational 
calculus on £7 gives the associated discrete Euler- Lagrange equation: 

Vkel,...,N-l, Q^- 2 Qk + Qk-i = _ W (Q fc ), Qe(R d f+\ 
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The proof of theorem [2] (|17J) is based on the following result: 

Lemma 1 (Discrete integration by parts). For any F and G in (R d ) N+1 , we have: 
N N-l 

^(A_F) fe G fe = F k(^+G) k + tUnQn - /oso). 



k=l 



k=0 



Discrete integration by parts emphasizes the asymmetric property of discrete operators 
A + and A_ which does not exist in the continuous space with the operator 4r. This 



characteristic leads to an asymmetry in ( EL^ ): indeed, we have a composition between 
the two discrete operators A + and A_. We notice that this asymmetry does not appear 
in the continuous space in (lELp . 



3.2 No coherence between direct discrete embedding and variational 
integrator 

We resume in figure [T] the previous processes in the setting of Gauss-FDEa. We are inter- 
ested in the coherence or not of this discrete embedding: do the direct discrete embedding 
and the variational integrator lead to the same numerical scheme for the Euler-Lagrange 
equation ? 



Lagrangian functional C ■ 



Functional discretization 



Discrete Lag. functional C? 



Variational principle 



■lational integral 



Euler-Lagrange equation ([EL]) 



Discr. var. principle 



Direct discrete embdedding 



Num. schemes (J22]) / $ELh§ 



Figure 1: No coherence between direct discrete embedding and variational integrator of the 
Euler-Lagrange equation in the setting of Gauss-FDEa. 



The previous study leads to a default of coherence of Gauss-FDEa. Indeed, algo- 
rithms obtained by direct discrete embedding (|2,2p and obtained by discrete variational 
principle do not coincide. A temporal asymmetry appears in ^EL^ ) which comes 

(mathematically) from lemma [TJ 
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The direct discrete embedding respects the law of semi-group of the differential operator 
A- 

df 

d d d direct / a \ / an a 9 

372 = 17 T77 ► (-aA a )o(-aA a ) = A;. 

dt CLt OLt discretization 

As it is well-known from the numerical analysis point of view, [9], the composition of 
the discrete operator — o"A CT provides unstable numerical schemes. We loose the order of 
approximation of differential operator by composition: indeed, A^ approaches J™ with an 
order 1 and not an order 2. 

On the contrary, a variational integrator is not based on the algebraic construction of 
the differential equation via the differential operator but mainly on a dynamical approach 
via the extremals of a functional. As a consequence, the discrete Euler-Lagrange equation 



( ELh ) does not necessary respect the algebraic property of composition. Nevertheless, the 
variational integrator uses the composition of discrete operators — A + and A_ and this 
composition provides an approximation of with an order 2. 



The problem is to understand why this asymmetry does not appear with direct dis- 
crete embedding ? It seems that we miss dynamical informations in the formulation of 
Lagrangian systems at the continuous level which are pointed up in the discrete space with 
the asymmetric discrete operators (— <tA ct ) ct=± . 

This default of coherence can be corrected using a different writing of the initial Euler- 
Lagrange equation. Indeed, we introduce in section H] differential operators which conserve 
the temporal asymmetry in the continuous space: namely the asymmetric derivatives. 



4 Discrete embedding of asymmetric Lagrangian systems 

The usual way to derive differential equations in Physics is to built a continuous model 
using discrete data. However, this process gives only an information in one direction of 
time. As a consequence, a discrete evaluation of the velocity corresponds in general at the 
continuous level to the evaluation of the right or left derivative. In general, we replace the 
right (or left) derivative by the classical derivative jL However, this procedure assumes 
that the underlying solution is differentiable. This assumption is not only related to the 
regularity of the solutions but also to the reversibility of the systems (the right and left 
derivatives are equal). In this section, we introduce asymmetric Lagrangian systems which 
are obtained with functionals depending only on left or only on right derivatives. We prove 
in this case that Gauss-FDEa is coherent. 
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4.1 Asymmetric notations and asymmetric integration by parts 

Like operators — A + and A_, we define differential operators which use only either future 
informations or past informations. They are just the minus right and the left classical 
derivatives: 

Definition 8. For f : [a, b] — > M. d smooth enough function, we denote: 

Vte[a,b[, d + f(t)= lim /^-f + k) 

and 

Vt€]o,6], d.f(t)= lim 

h->0+ h 

Of course, for a differentiable function /, we have d-f = —d + f = f. However, as we are 
going to see in the next subsection, it is interesting to use these notations in order to keep 
dynamical informations. 

Lemma 2 (Asymmetric integration by parts). For f,g : [a,b] — > M d smooth enough 
functions, we have: 

1 d-f(t)g(t)dt = f f(t)d + g(t)dt + f(b)g(b) - f(a)g(a). 



Let us notice that the asymmetric integration by parts is the continuous analogous of the 
discrete integration by parts (see lemma [1]). 

4.2 Asymmetric Lagrangian systems 

In this section, we introduce asymmetric Lagrangian systems which take into account the 
temporal asymmetry by using the differential operator —ad a . 

Definition 9. An asymmetric Lagrangian functional is an application: 

C 7 : C 2 ([a,b],R d ) — > M 

q i — ► J^L(q(t),-o-d a q(t),t)dt 



where L is a Lagrangian. 

Then, we obtain the following characterization of the extremals of an asymmetric La- 
grangian functional as solution of an asymmetric differential equation: 

Theorem 3 (Variational principle). Let C a be an asymmetric Lagrangian functional as- 
sociated to the Lagrangian L and let q € C 2 ([a,b],M. d ). Then, q is an extremal of C 7 if and 
only if q is solution of the asymmetric Euler- Lagrange equation: 

Vt€]a,6[, §(q(t),-ad a q(t),t)-ad. a (^(q,-ad a q,.))(t)=0. (EL a ) 
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Indeed, by using the asymmetric integration by parts in the proof of theorem [3l the asym- 
metry appears in ( EL a ). Its origin is the continuous analogous of the asymmetry in ( ELh ). 
Now, we are interested in discrete embeddings of the asymmetric Euler-Lagrange equation. 



4.3 Direct discrete embeddings and variational integrators of asymmet- 
ric Lagrangian systems 



In order to discretize flEXg-D , we have to discretize two differential operators at the same 
time. In the discrete space, we replace the differential operator d+ (respectively d_) by 
A + (respectively by A„). Then, we obtain an asymmetric version of the Gauss finite 
differences embedding mixing A + and A_: 

Definition 10. We call the asymmetric Gauss-FDEa embedding the definition of the 
following elements: the application 

disc: C°{[a,b],R d ) — > (R d ) N+1 

q i — > (q(t k ))k=o,...,N 
the a -quadrature formula of Gauss and the discrete operators 

A+ : (R d ) N+1 f^ N 

I y ( Qk—Qk + 1 \ 

k=0,..,N-l 



r\ _ i r\ \ . ( Qk-Qk+i 

H — Wk)k=0,...,N 1 > I h 



and 

A_: (R d ) N+1 — > (if 

Q = (Qk)k=o,...,N 1 — > (°*~?*~ 1 ) fc=1 N 

These discrete operators are respectively the discrete versions of the operators d+ and d- 
in the asymmetric Gauss-FDEa. 

Definition 11. Let L be a Lagrangian and \EL a ) the asymmetric Euler-Lagrange equation 
associated. The direct discrete embedding of <\EL a ) with respect to the asymmetric Gauss- 
FDEa is given by: 

— (Q, -aA a Q, r) - aA_ a ( — (Q, -aA a Q, r) = 0, Q € (M d f+\ 

Lat us notice that the direct discrete embedding of the asymmetric Euler-Lagrange equation 
gives the discrete Euler-Lagrange equation ( ELh\) (see definition [2]). Obviously, as the 
asymmetric Gauss-FDEa gives the same discrete Lagrangian functional as the symmetric 
case, we deduce the coherence of this discrete embedding. 

We notice that the asymmetric Gauss-FDEa unifies the algebraic approach and the 
dynamical approach in order to discretize Lagrangian systems. Indeed, the discrete and 
continuous versions of the asymmetric differential operators satisfy the same algebraic 
properties (temporal asymmetry, asymmetric integration by parts,...). 
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5 Introduction to fractional calculus 



5.1 Notions of Griinwald-Letnikov and Riemann-Liouville 

Fractional calculus is the generalization of the derivative notion to real orders. There are 
many ways generalize this notion. We can generalize the derivative of a monomial and 
then generalize the derivative of functions which can be written in power series. We can do 
the same thing by generalizing derivative of an exponential or again derivative of cosinus 
and sinus. We can also use a generalization of the derivative by using a convolution. For 
more details, we refer to [19J. 

In this paper, we are going to use the classical notions of Griinwald-Letnikov and 
Riemann-Liouville. The following definitions and results are extracted from 



• Fractional derivatives of Griinwald-Letnikov The notion of Griinwald-Letnikov 
comes from this simple result proved by induction: 

let / be a smooth enough function defined on [a, b] and let n £ N, then: 



where CZ 



1 n 

Vte}a,b], fW(t)=lim + -^2(-iyc:f(t-rh) 
n{n — l)...(n — r + 1) 



h->o+ h n 

r=0 



From this formula, Griinwald-Letnikov obtains this following generalization: 

Definition 12. Let a > and let f be an element o/C n+1 ([a, b], W 1 ) where n = [a]. 
Griinwald-Letnikov fractional left derivative of order a with inferior limit a of 

f is: 

n 

Vt€]a,b], D a _f(t)= lim -J2(-l) r C r J(t-rh). 

nh=t-a r=0 

Griinwald-Letnikov fractional right derivative of order a with superior limit b 

of f is: 

1 n 

Vt€[a,6[, D%f(t)= lim + rh). 

nh=b-t r=0 

• Fractional derivatives of Riemann-Liouville The notion of Riemann-Liouville 

comes from this other simple result proved by induction: 

let / be a continuous function on [a, b], let t in [a, b] and let n G N, then: 

/ " f( tn )dt n ...dt 1 = T ^— f (t-yr- 1 f(y)dy. 
{n - 1)! J a 



a J a J a 



From this formula, Riemann-Liouville obtains a generalization of the integral notion: 
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Definition 13. Let a > and let f be a continuous function. 

Riemann-Liouville fractional left integral of order a with inferior limit a is: 

Vie [a, 6], RL DZ a f(t):=-l- [ \t - y)^ 1 f(y)dy. 

r(a) J a 

Riemann-Liouville fractional right integral of order a with superior limit b is: 

Vt € [a, b], RL D- a f{t) : = -L [\y - tr~ l f(y)dy. 

where T is the Gamma function which is an extension of the factorial function. 

Finally, Riemann and Liouville obtain a generalization of the derivative notion by derivating 
fractional integrals: 

Definition 14. Let a > and let f be an element o/C n+1 ([a, b], M. d ) where n = [a]. 
Riemann-Liouville fractional left derivative of order a with inferior limit a is: 

jn+l , . 

Vt e]a,b], RL D«f(t) := 1^ ( RL DT in+1) f) (t). 

Riemann-Liouville fractional right derivative of order a with superior limit b 

is: 

Vt G [a, b[, RL D°f(t) := ^ ( ^ X) /) (*)■ 
5.2 Links between these different notions 

Although the two approaches are different, they are linked by the following result, [22j: 

Theorem 4. Let a > 0, a < b two reals and n = [a]. Let f be a function in C n+l ([a, b], M. d ) 
then the notion of fractional derivative of Riemann-Liouville and the one of Griinwald- 
Letnikov coincide. So we denote: 

RL D«f(t) = Dtf(t) 

and 

RL D a + f(t) = D%f(t). 

Now, the following result proves that the Riemann-Liouville 's and the Griinwald-Letnikov's 
fractional derivative notion is actually an extension of the classical derivative notion: 

Proposition 1. Let n € N* and f be an element of C n+1 ([a, b], M ). Then, the following 
equalities hold: 

Dlf = and D\f = (-1)"/H 

There exist fractional left derivatives and fractional right derivatives: this notion takes into 
account the temporal asymmetry. Moreover, there exists a formula of fractional integration 
by parts: 
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Lemma 3 (Fractional integration by parts). Let < a < 1 and let f and g be two functions 
in C 1 ([a, b], We suppose that f(a) = f{b) = or g(a) = g(b) = 0. Then, we have: 

f D°Lf{y)g{y)dy= f ' f{y)D%g{y)dy. 

J a J a 

Let us notice that the fractional integration by parts is the fractional analogous of lemma 
[2 it transforms right fractional derivative into left fractional derivative. We will see that 
we have same behaviors in the discrete space with the discretized fractional operators (see 
lemma [4| . 



6 Discrete embeddings of fractional Lagrangian systems 

In recent years, an important activity has been devoted to fractional Lagrangian systems 
for the purpose of optimal control, mechanics, engineering and Physics ([2], [3], [16]). There 
also exist many studies concerning the discretization of fractional differential equations ([2], 

®, ®, m, 0). 

In this section, following our previous approach on classical Lagrangian systems, we 
define fractional discrete embeddings and fractional variational integrators. We prove that 
a fractional discrete embedding is coherent. 

6.1 Fractional Lagrangian systems 

Definition 15. A fractional Lagrangian functional of order < a < 1 is an applica- 
tion defined by: 

£ a : C 2 ([a,b],m d ) — > R (6.1) 
q ,— > f*L(q(t),-o-D%q(t),t)dt 

where L is a Lagrangian. 

We can give a characterization of extremals of a fractional Lagrangian functional as solu- 
tions of a fractional differential equation: 

Theorem 5 (Variational principle). Let C 7 be a fractional Lagrangian functional of order 
< a < 1 associated to the Lagrangian L and let q be an element of C 2 ([a, b], W 1 ). Then, 
q is an extremal of C a if and only if q is solution of the fractional Euler- Lagrange 
equation: 

Vi G]a,&[, §(q(t), -aD%q(t),t) - aD% (§(q, -aD^q, .)) (i) = 0. (EL f ) 
As in the asymmetric case, by using the fractional integration by parts, we obtain a frac- 



tional differential equation which presents a temporal asymmetry. Indeed, in {ELf), we 
have a composition of the D°L operator and the D°_ operator. 
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6.2 Direct Gauss Griinwald-Letnikov embeddings of fractional Lagrangian 

systems 

We are interested in discrete embeddings of fractional Lagrangian systems. By referring to 
the notion of Griinwald-Letnikov (|15j), we give the following definition: 

Definition 16. The Gauss Griinwald-Letnikov embedding denoted by Gauss-GLEa 
is the definition of the following elements: the application 

disc: C°([a,b],R d ) — >• (R d ) N+1 

q i — > (q{t k )) k =o,...,N 
the a -quadrature formula of Gauss and the discrete operators 

A a . (R d ) N + 1 _> (R d ) N 



and 



Q = (Qk)k=0,...,N 



\ r=0 



f I N - k r \ 

Q = (Qk)k=0,...,N 1 > ( ^ (- l ) rC aQk+r ] 

V r=0 J k=0,..,N-l 

These discrete operators are respectively discretized operators of the operators D°_ and D" . 

With such discrete operators, we have a formula of discrete fractional integration by parts 
which is exactly the discrete analogous of the fractional integration by parts at the contin- 
uous level (see lemma [3j): 

Lemma 4 (Discrete fractional integration by parts). Let F and G be elements of {M. d ) N+1 
such that Fq = Fn = or Gq = Gn = 0. Then, we have: 

N N-l 
fe=l fc=0 

Proof. The proof of this lemma is only based on the inversions of the two sums and on 
substitution of variables in the sums. We prove the lemma only in the case Fq = Fjf = 0. 
The case Go = Gn = is similar. 

As .F/v = 0, we can sum k from to N and we obtain by definition of (A"G)k'- 

J2k=0 Fk{^%G) k = ^J2k=oT,r=O k (- 1 Y C 'a F kGk+r 

1 v^-ZV v-*iV— 77 -\\rr<rT?n 

— W 2^r=0 2^k=0 \ L ) ^a p k^k+T 

_ 1 v^./V v^Af (_-\\r/^rTp ri 

— h a <L>r=0 l^k=r\ L ) ^a r k-r^k 
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Then, as Fq = 0, we can sum k only from 1 to N and we obtain: 



□ 



As in the asymmetric case, we notice that fractional operators and their discretized take 
into account the temporal asymmetry. The integration by parts on these operators leads 
to the same behavior: the fractional right derivative is transformed into the fractional left 
derivative. Then, we can expect the coherence of Gauss-GLEa. Now, let us give the 



definition of the direct discrete embedding of (ELf) with respect to Gauss-GLEa: 



Definition 17. Let < a < 1 and let L be a Lagrangian and (ELf ) the fractional Euler- 
Lagrange equation associated to L. The direct discrete embedding of (ELf) with respect to 
Gauss-GLEa is given by: 

—(Q,-aA«Q,T)-aA%(—(Q,-aAZQ,T))=0, Q G (R d ) N+1 . (6.2) 

6.3 Discrete calculus of variations on discrete fractional Lagrangian func- 
tional 

Giving Gauss-GLEa allows us to define the discretization of fractional Lagrangian func- 
tionals: 

Definition 18. Let < a < 1 and let C a be the fractional Lagrangian functional associated 
to the Lagrangian L. The discrete fractional Lagrangian functional associated to C a 
with respect to the embedding Gauss-GLEa is given by: 

C u h : (R d ) N+1 — > R 

Q = (Qk)k=o,...,N 1 — > h L(Q k , (-aA"Q) k ,t k ), 

where I + = {0, ...,N — 1} and J_ ={1,...,N}. 

Discrete extremals of can be characterized as solution of a system of equations: 

Theorem 6 (Discrete variational principle). Let < a < 1 and let be a discrete 
fractional Lagrangian functional associated to the Lagrangian L with respect to the discrete 
embedding Gauss-GLEa. Then, Q in (IR^)^"^ 1 is a discrete extremal of if and only if 
Q is solution of the following system of equations, called the discrete fractional Euler- 
Lagrange equation: 

r)T BT 

— (Q,-aA a a Q,T) - aA a _ a (—(Q,-aA«Q,T)) = 0, Q G (R d ) N+1 - (EL fth ) 
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Proof. We write the proof only in the case a = —. The proof in the case a = + is similar. 
Let Q in (M d ) Ar+1 and W a discrete variation of (JR d ) 7V+1 . Then: 

N 

C h (Q + W) = hY,L(Qk + W k ,(A a _Q) k + (A a _W) k ,t k ) 

k=l 

- BL 81 
= Cj;(Q) + h^2— (Q k ,(A a _Q) k ,t k )W k + —(Q k ,(A«Q) k ,t k )(A a _W) k 

k=l X V 

+o(\\W\\) 

where o is the Landau's notation. Thus, by definition of the differential, we obtain: 

V dL,„ dL 



DC h {Q){W) =hJ2 -Q-(Qk, (A a _Q) k ,t k )W k + — (Q k , (A a _Q) k ,t k )(A«W) k 

k=l X V 

Since W is a discrete variation of (M d ) Ar+1 , Wq = Wn = 0. By applying the discrete 
fractional integration by parts on the second term of the previous sum, we obtain: 

DCl(Q)(W) = 0^h^( — (Q k , (A a _Q) k ,t k ) + (A° (— (Q, A°LQ, r))) k )W k = 0. 

k=l ^ X V ' 

Thus, Q is a discrete extremal of if and only if Q is solution of the discrete fractional 



Euler-Lagrange equation (ELfh) as defined in theorem El 

□ 

As we expected before, we notice that we have coherence of the discrete embedding Gauss- 



GLEo of the fractional Euler-Lagrange equation. Indeed, formulas (|6.2p and ( EL t ^ ) 
coincide. 



Conclusion 

In this paper, we have introduced the notion of discrete embeddings and constructed the 
corresponding variational integrators in the classical and fractional cases. The discrete 
embedding has to conserve the integration by parts on the discrete level in order to be 
coherent. 

The Gauss finite differences embedding or Griinwald-Letnikov embedding are only here 
as illustration of our point of view. In order to construct more efficient numerical schemes, 
one can choose appropriate discretization for the differential operators and also more elab- 
orate quadrature formulas for integrals. In particular, we can have a look in |15] for special 
discretization of the fractional derivatives. 



In |23| . Riewe has introduced fractional Lagrangian functionals in order to model dis- 
sipative effects. However, his formalism does not provide an equivalence between solutions 
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of a dissipative equation and critical points of fractional Lagrangian functional. Moreover, 
his construction is based on very stringent assumptions on functions. The asymmetric 
fractional calculus of variations (where the functional space is split in two) introduced in 
|13| solves this problem and provides a full equivalence between critical points of an asym- 
metric fractional Lagrangian functional and solutions of a dissipative system (as the linear 
friction case or the diffusion equation for example). A natural extension of our work is to 
extend the fractional discrete embedding to cover the asymmetric fractional case. 

Several difficulties arise in the numerical study of the convection-diffusion equation, 
in particular the presence of non physical numerical oscillations for high values of the 
Reynolds number. An asymmetric fractional Lagrangian formulation of the convection- 
diffusion equation is obtained in [12| . Our idea is to construct an adapted fractional 
variational integrator in order to control the discrete solutions of the numerical scheme. 
This will be done in a forthcoming paper. 
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